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ABSTRACT 

The rotation velocity asymmetry v observed in spiral galaxy HI rotation curves linearly cor- 
relates with the effective potential force from the 10 closest neighboring galaxies normalized for 
the test particle mass and the gravitational constant. The magnitude of the potential force from 
a close galaxy is proportional to the luminosity of the close galaxy and inversely proportional to 
the square of the distance from the close galaxy to the target galaxy. The correlation coefficient 
is 0.99 and F test is 0.99. Also, the slope of the rotation curve in the disk region of a galaxy 
from rising to flat to declining is qualitatively correlated with increasing asymmetry and, hence, 
to the net force from other galaxies. The result is based on a sample of nine spiral galaxies with 
published Cepheid distances and rotation curves and with a wide range of characteristics. These 
relationships are interesting not only for their predictive power but also because (1) they suggest 
a galaxy's dynamics and the shape of its rotation curve are related to the potential energy and 
distance of neighboring galaxies, (2) they suggest a rising rotation curve in the disk region is 
intrinsic, (3) they suggest the Tully-Fisher relationship's assumption that the mass to intrinsic 
luminosity ratio is constant among galaxies is valid, and (4) they are inconsistent with MOND, 
dark matter, and the Linear Potential Model. 

Subject headings: cosmology :theory — galaxies:fundamental parameters — galaxies:kinematics and dy- 
namics 



1. INTRODUCTION 

Asymmetry is often observed in rotation curves 
of spiral galaxies. Like the rotation curve, asym- 
metries indicate the forces and dynamics in a 
galaxy. Small deviations in the rotation curves of 
a few kms~ l are well known (Shane and Bieger- 
Smith 1966). Large deviations are less appreci- 
ated (Jog 2002) because the observational data are 
generally averaged. Hence, only highly asymmet- 
ric cases such as NGC 4321 (Knapen et al. 1993) 
and NGC 3031 (Rots 1975) are recognized. The 
rotation curves of all the nearby galaxies where 
the kinematics are studied have asymmetry such 
as NGC 0224 (Simien et al. 1978) and NGC 0598 
(Colin and Athanassoula 1981). Asymmetry has 
also been reported in the inner parts of the opti- 
cal disk (Sofue 1988). Rotation curve asymmetry 
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appears to be the norm rather than the excep- 
tion (Jog 2002). Weinberg (1995)and Jog (1997) 
proposed the implied mass asymmetry is due to 
a galaxy interaction. In this Paper, the rotation 
curve asymmetry from our viewpoint is related to 
the strength of the forces from neighboring galax- 
ies. A qualitative correlation of rotation curve 
shape with rotation curve asymmetry is shown. 
This result suggests the rotation curve shape is 
naturally rising and neighboring galaxies cause flat 
and declining rotation curves. These results are 
inconsistent with MOND, dark matter, and the 
Linear Potential Model. 

2. MODEL 

Hodge & Castelaz (2003) (hereinafter Paper 1) 
found the observed magnitude of a target galaxy 
depends on the radiant energy of neighboring 
galaxies and on the intrinsic luminosity of the tar- 
get galaxy. The potential energy field of neighbor- 
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ing galaxies may also produce observable effects in 
target galaxies. 

Consider a particle in the disk region of a spiral 
galaxy. If there are no outside forces, the intrinsic 
orbit of the particle will be achieved by a delicate 
balance of the intrinsic forces of the galaxy. The 
circular orbits of galaxy disks (Binney & Merri- 
field (1998), pages 723-4) imply the intrinsic forces 
are centrally directed. Therefore, a small, external 
force directed nearly uniformly across the plane of 
the galaxy may cause a perturbation in the orbit 
and in the rotational velocity of particles in the 
disk. Thus, the asymmetry of the rotation velocity 
curve may be a sensitive indicator of the external 
forces exerted on a galaxy's particles. 

To test the effects of the potential energy field of 
neighboring galaxies, a measure of galaxy mass is 
needed. Commonly, mass is derived from rotation 
curve models. However, such mass calculations are 
lacking for most galaxies. Paper 1 and Aaronson, 
Huchra, & Mould (1979) suggested the success of 
the Tully-Fisher relationship (TF) implies that the 
ratio of intrinsic mass to intrinsic luminosity is 
constant among galaxies. Therefore, 

mass cx lO"^/ 2 - 5 , (1) 

where M b /j\ (in units of magnitude) is the intrinsic 
absolute B band magnitude M b of the j th neighbor 
galaxy. The adjective "intrinsic" means after the 
luminosity effects of neighbor galaxies have been 
subtracted (Paper 1). 

The force on a particle in the k th target 
galaxy exerted by a number n of neighboring close 
galaxies is, 

n rt 

F (fc) =£G m mlO- M ^/ 2 - 5 9M, (2) 

i=i n m 

where Ruk) 1S the vector distance from the k th 
target galaxy to the j th neighbor galaxy, G m is 
the constant of proportionality such that the units 
of F( fe ) are dynes, and m is the mass of the test 
particle in the target galaxy. The bold letters de- 
note vectors. The same letter not bold denotes 
the magnitude of the vector. The letters in paren- 
theses are indices denoting the galaxy. For ease of 
calculation M b ^ is without the correction found 
in Paper 1. Hence, the G m factor also contains a 
constant which is the average change in luminosity 



of the j galaxy due to re-emitted radiant energy 
from neighboring galaxies. From Paper 1, the de- 
viation from average of the re-emitted energy will 
produce an error in M b (j\ of ± 0.3 magnitude. 

Define F m ^ as F( k ) in a coordinate system 
that is galactocentric, right-handed, and Carte- 
sian, 

F m (k) — {.Fminor(k) ~t~ F ma j or (k) Fp [ ar (k)) i 

(3) 

where F minor ^ k ) , F major ( k ) , and F polar ( k) are the 
component vectors of Ft^ in the direction of the 
minor axis, major axis, and polar axis of the k th 
target galaxy, respectively. Since the F ma j or ( k ^ 
acts radially along the major axis, the effect of the 
central force of a galaxy F c along the major axis 
is reduced if the measurement of asymmetry is at 
the same radius from the center on each side of the 
galaxy. Define the rotation curve asymmetry Vd as 
the maximum difference in rotation velocity along 
the major axis between one side of the center of a 
galaxy to the other measured at the same radius 
from the center of the galaxy. 

Consider the effect of F m ^) on Vd is exerted 
proportionally by a vector K, 

v d = \K • F m{k) \, (4) 

= K rn i nor e rn i nor ^K rna j or G rna j or -\- Kpolar&polar ■> 

(5) 

where e m i nor , e ma j or , and epoiar are the unit vec- 
tors in the minor axis direction, the major axis 
direction, and the polar axis direction of the k th 
target galaxy, respectively, K mmor , K major , and 
Kpoiar are coefficients of the corresponding unit 
vectors, and the | | means the absolute value. 
The direction of e ma j or is directed eastward along 
the axis defined by the position angle (PA). The 
e m j nor is perpendicular to e ma j or , in the plane 
of the galaxy, and directed away from the Milky 
Way. A right-handed coordinate system means 
e m inor x emajor = e po iar- In conformity to Paper 
1, a positive orientation sign implies e minor has 
a component directed southward and e po i ar has 
a component directed away from the Milky Way. 
A negative orientation sign implies e minor has a 
component directed northward and e po i ar has a 
component directed toward the Milky Way. 

The kinematic minor axis is the direction along 
which the line-of-sight motion equals the systemic 
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velocity. The kinematic major axis is the direc- 
tion along which the line-of-sight motion peaks as 
a function of angle about the galactic center. The 
tilted ring model predicts the kinematic minor axis 
is perpendicular to kinematic major axis. Galax- 
ies with warped disks may have curved kinematic 
axes. In this Paper, the major and minor axes are 
not kinematic axes when a warp is present. 

The K is found by plotting Vd versus \K»F m ^ | 
such that a best fit of equation (4) results. 

Since the distance to neighboring galaxies is 
much greater than the diameter of galaxies, the 
first approximation is that F m r k \ is uniform across 
a galaxy. The asymmetry in the HI rotation curve 
is a sensitive tracer of \K • F m ^\. 

3. DATA AND ANALYSIS 

The criteria for choosing target galaxies for 
the analysis are (1) the galaxy has a published 
Cepheid distance D c and a published HI rotation 
curve asymmetry, (2) the distance to the galaxy 
must be large enough that the Milky Way's contri- 
bution is negligible, (3) the distance to the galaxy 
must be large enough that the contribution of pe- 
culiar velocity to the redshift z measurement is rel- 
atively small among the galaxies around the target 
galaxy, and (4) a published apparent magnitude in 
the B band to& must be available for both the tar- 
get and close galaxies. Nine galaxies were found 
to meet these criteria. 

The data for the sample of nine target galax- 
ies are presented in Table 1. The morphology 
was taken from the NED database 3 . The rota- 
tion curve shape and measured Vd were taken from 
the references. The galaxy orientation sign which 
was used to determine the component vectors of 
F m (k) was taken from Paper 1. The total apparent 
corrected B-band magnitude m& is "btc" in the 
LEDA database 4 . The value of "btc" for close 
galaxies was used since I band data is generally 
unavailable. The D c was taken from Freedman et 
al. (2001) except as noted in Table 1. 

The sample of target galaxies has low sur- 
face brightness (LSB), medium surface brightness 

3 The NED database is available at: 

http://ncdwww.ipac.caltcch.edu. 

4 The LEDA database is available at: http://leda.univ- 

lyonl.fr. 



(MSB), and high surface brightness (HSB) galax- 
ies, galaxies with a range of 21-cm line width at 
20% of peak value corrected for inclination Wfo 
from 122 fcms -1 to 654 kms^ 1 , includes LINER, 
Sy, HII and less active galaxies, galaxies with a 
absolute B magnitude ("mabs" from the LEDA 
database) from -22.11 mag. to -17.03 mag., a dis- 
tance range from 1.13 Mpc to 14.73 Mpc, and field 
and cluster galaxies. 

The NED database was used to assemble a list 
of galaxies closest to each of the target galax- 
ies. The selection of close galaxies used the 
Hubble Law with a Hubble constant H a of 70 
kms~ 1 Mpc~ 1 and galactic z to calculate the dis- 
tances to galaxies. The angular coordinates were 
taken from the NED database. For each target 
galaxy, the distance from each target galaxy to 
each of the other galaxies was calculated. The 
neighbor galaxies with the smallest distance were 
chosen as the close galaxies. 

The uncertainty of the value of H a is high. The 
use of H Q herein is restricted to only choosing the 
close galaxies to each target galaxy. As the relative 
distances become larger, the distance uncertainty 
increases if only because of a greater variation in 
peculiar velocity. Therefore, there is a practical 
upper limit on the number of galaxies that can be 
used to evaluate F m ^ k y If too few close galaxies 
are used, F m r k \ will be underestimated. Following 
Paper 1, the number of close galaxies found to 
produce the highest correlation coefficient was 10. 

The rotation curve of NGC 4321 at 4 arcmin 
implies it has a declining shape. NGC 4321 has 
a large asymmetry in the south west lobe of 209 
fcms -1 which cannot be interpreted as disk rota- 
tion (Knapen et al. 1993). The Ha rotation curve 
is flat out to 20 kpc with Vd — 52Akms~ 1 (Arse- 
nault et al. 1988). The receding side of the HI ro- 
tation curve of NGC 4321 is declining and appears 
more normal. The approaching side of the HI ro- 
tation curve of NGC 4321 is "lopsided" (Knapen 
et al. 1993) and is sharply rising after a slight de- 
cline. Guhathakurta et al. (1988) suggested the 
HI rotation curve is declining. 

The R(jk) and M b ^ in equation (2) were calcu- 
lated as in Paper 1. Because M b tj) is more uncer- 
tain without a secure distance (D c ) to each galaxy, 
the value used for M b ^ was not the intrinsic mag- 
nitude. Therefore, the value of F m ( k - ) was overes- 
timated. This causes the value of \K • F m ( k )\ to 
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have a larger scatter. 

Figure 1 is a plot of v d versus \K • F m ^)\- The 
coefficients of K were determined such that vari- 
ance between the line with a slope of one and data 
(F test) is maximized and the sum of the difference 
between the points and the line is minimized, 

v d = S\K.F m{k) \ +1, (6) 

where 

K = (9.70 x l(T 7 )e mmor + (1.88 x \Q- 7 )e major 
-(4.24 x l(T 8 )e po;al . 

fcms -1 dyne~ 1 G m m, (7) 

K = 9.89 x 10~ 7 km s~ x dyne~ 1 G m m, the slope 
S = 1.00±0.06, and the intercept / = Q±'3km S - 1 . 
The error is 1 er, the correlation coefficient is 0.99, 
and the F test is 0.99. 

Figure 2 is a plot of Vd versus rotation curve 
type. Figure 2 suggests a qualitative correlation 
between w d and the slope of the sample galaxies' 
rotation curve in the disk region. 

4. DISCUSSION 

The v d is defined at the radius from the center 
along the major axis with the maximum asymme- 
try in rotation velocity. We speculate this implies 
the dynamics of galaxies are a function of radius 
from the center rather than a number of scaling 
factors. 

Since both Paper 1 and this Paper require the 
orientation sign of the target galaxy to be estab- 
lished, the method of this Paper may confirm the 
orientation sign that the Paper 1 method finds. 

Paper 1 used luminosity (Mb) as an indicator of 
the radiant energy from a neighbor galaxy. This 
Paper used luminosity as an indicator of mass in 
a neighbor galaxy. Luminosity, mass, and radi- 
ant energy of a galaxy appear tightly related. The 
difference between radiant energy's action and po- 
tential energy's action on a target galaxy is in the 
way energy is transferred. The amount of radiant 
energy transferred to a target galaxy depends on 
the surface area presented to the neighbor galaxy. 
Hence, the radiant energy transferred is maximal 
if the neighbor galaxy is along the target galaxy's 
polar axis and minimal if the neighbor galaxy is 
in the plane of the target galaxy. The potential 



energy transferred is at a maximum if the neigh- 
bor galaxy is in the target galaxy's plane and at a 
minimum if the neighbor galaxy is along the polar 
axis of the target galaxy. 

Paper 1 concluded that the contribution of the 
conservative, potential field of neighboring galax- 
ies contributed no energy to the B-band luminosity 
because matter in a galaxy orbited in radial shells 
of constant mass. The target galaxy's movement 
in the intergalactic potential field could cause only 
minimal heat change. The correlation coefficient 
between the Paper 1 correction factor C/ and v d is 
-0.04, between C/ and \K»F m ^) \ is -0.09, between 
the Paper 1 impinging energy E a and v d is 0.61, 
and between E a and \K • F m r k \\ is 0.65. These 
poor correlations imply the potential energy field 
of the close galaxies has little or no effect on the 
Mb measurement of the target galaxies. 

Persic, Salucci, & Stel (1996); Persic & Salucci 
(1990); Sofue et al. (1999) found a strong corre- 
lation between a galaxy's luminosity and rotation 
curve shape. Low luminosity galaxies have rising 
rotation curves; high luminosity galaxies have flat 
and declining rotation curves. This Paper found 
close neighbor galaxies appear to cause higher 
asymmetry in the rotation curve and to cause flat 
and declining rotation curves. Also, this paper 
suggests that without the influence of neighbor- 
ing galaxies (intrinsic), the rotation curve would 
be rising. Paper 1 found close neighbor galaxies 
appear to cause higher luminosity in a galaxy de- 
pending on the orientation of the galaxy relative to 
the close galaxies. Therefore, neighboring galax- 
ies cause differences in the mass/luminosity ratio 
among galaxies by increasing the apparent lumi- 
nosity as outlined in Paper 1. 

The dark matter model (DMM) expects the in- 
trinsic rotation curve to be declining. If dark mat- 
ter (DM) is considered part of a galaxy's mass, the 
constant intrinsic mass to luminosity ratio among 
galaxies implies the ratio of DM to luminous mass 
is also constant among galaxies. To have a rising 
rotation curve caused by DM implies the poten- 
tial caused by the DM must increase with radius 
while the luminous mass is decreasing with radius. 
If galaxies have differing mass because of differing 
mass infall rates, then the constant ratio of DM to 
luminous mass and the required DM distribution 
seems implausible. 

The linear potential model (LPM) explored by 
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Mannheim and Kazanas (Mannheim & Kmetko 
1996) proposes a gravitational potential of the 
form V(r) = -f3/r + jr/2. The LPM predicts 
rising rotation curves are the intrinsic shape and 
flat rotation curves are a result of mass distribu- 
tion within the galaxy. Also, the LPT predicts 
declining rotation curves require other galaxy in- 
fluence with a linear potential and a very deli- 
cate balance. The modified Newtonian dynamics 
(MOND) model (Bottema et al. 2002) expects the 
intrinsic rotation curve to be flat and flat rotation 
curves to be axis symmetric. The finding of this 
Paper is inconsistent with the LPM, the MOND 
model, and the DMM. 

5. CONCLUSION 

The rotation curve asymmetry is linearly 
related to the force from neighboring galax- 
ies. The relation is v d = 1.00(±0.06) \K • 
F m (k)\ + 0(±3) at 1 a where K = (9.70 x 
lCT 7 )e minor + (1.88 x 10- 7 )e major - (4.24 x 
10~ 8 )e po / ar kms" 1 dyne~ 1 G m m. The correlation 
coefficient is 0.99 with F test of 0.99. 

Rising rotation curves have little rotation curve 
asymmetry and flat and declining rotation curves 
have greater asymmetry. 

Paper 1 is correct in assuming the potential 
field of neighboring galaxies has no measurable im- 
pact on the Mt, of target galaxies. 

The orientation sign adopted for Paper 1 is con- 
firmed for the galaxies in this Paper. This includes 
the problematic galaxies in Paper 1 NGC 3319 and 
NGC 4535. 

This research has made use of the NASA/IPAC 
Extragalactic Database (NED), which is operated 
by the Jet Propulsion Laboratory, California In- 
stitute of Technology, under contract with the Na- 
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|K • F | (kms 1 ) 

m (It) 

Fig. 1. — Plot of Vd (kms^ 1 ) versus \K • 
Fm(k)\ (kms^ 1 ). The equation of the line 
is Vd = \K • F m (fc)| where K = (9.70 x 
lCT 7 )e minor + (1.88 x I0- 7 )e major - (4.24 x 
10~ 8 )e po / ar km s^ 1 dyne~ 1 G m m. The line is a 
best fit of the data. 
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Rising Flat Declining 
RC type 

Fig. 2.— Plot of Vd (km s^ 1 ) versus rotation curve 
type. NGC 4321 has v d = 117fcms _1 . It's shape 
is either flat (Arsenault et al. 1988) or declining 
(Guhathakurta et al. 1988). 
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Table 1 

Data for the chosen galaxies. 



Galaxy 


NED typo a 


sign b 




RC d 


v d 6 


rof f 


\K.F m(k) \ * 


NGC 0925 


SAB(s)d;HII 




9.13±0.17 


D 


52±2 


1 


48±20 


NGC 2403 


SAB(s)cd 




3.14±0.36 


F 


8±2 


2 


14±10 


NGC 2841 


SA( r)b;LINER Sy 


+ 


14.07±1.57 h 


D 


16±2 


3 


25±10 


NGC 3031 


SA(s)ab;LINER Syl.8 




3.55±0.09 


D 


18±4 


4 


8± 4 


NGC 3109 


SB(s)m 


+ 


1.13±0.12' 


R 


1±2 


5 


2± 2 


NGC 3198 


SB(rs)c 




13.69±0.50 


F 


6±3 


6 


6± 3 


NGC 3319 


SB(rs)cd;HII 




13.44±0.57 


R 


6±2 


7 


15± 8 


NGC 4321 


SAB(s)bc;LINER HII 




14.33±0.46 


D 


117±2 


8 


121±60 


NGC 4535 


SAB(s)c 




14.80±0.32 


D 


50±2 


9 


49±20 



a Galaxy morphological type from the NED database. 

b The sign is the orientation of the polar axis of the target galaxy. A "+" sign indicates the polar axis is directed away from the 
Milky Way. The sign is taken from Paper 1. 

c Distance in Mpc from Ccphcid data from Frccdman ct al. (2001), unless otherwise noted. 

d Galaxy's HI rotation curve type according to slope in the outer disk region. R is rising, F is flat, and D is declining. 

e The maximum difference in rotation velocity (fcms -1 ) in the disk region along the major axis from one side of the galaxy to the 
other at the same radius. 

f References are for the rotation curve data for each galaxy. 

s In units of /cms -1 . 

h Distance is from Macri et al. (2001). 

distance is from Fcrrarcsc ct al. (2000). 

References.— (1) Krumm & Salpeter (1979); (2) Shostak (1973); (3) Fillmore, Boroson, & Dressier (1986); (4) Gottesman, Wcliachew 
(1975); (5) Jobin & Carignan (1990); (6) van Albada ct al. (1985); (7) Moore & Gottesman (1998); (8) Knapcn ct al. (1993); (9) 
Chincarini & de Souza (1985) 
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